Suppression of critical temperature in Ba(Fei_. c Co. r )2As2 with point defects introduced 

by proton irradiation 
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We report the effect of 3 MeV proton irradiation on the suppression of the critical temperature 
T c in Ba(Fei_ a! Co I! )2As2 single crystals at under-, optimal-, and over-doping levels. We find that 
T c decreases and residual resistivity increases monotonically with increasing dose. We also find no 
upturn in low-temperature resistivity in contrast with the a-particle irradiated NdFeAs(0,F), which 
suggests that defects induced by the proton irradiation behave as nonmagnetic scattering centers. 
The critical scattering rate for all samples estimated by three different ways is much higher than 
that expected in s±-pairing scenario based on inter-band scattering due to antiferro-magnetic spin 
fluctuation. 

PACS numbers: 74.62. En, 74.25.fc, 74.70.Xa 



Since the discovery of the high-T c iron-based super- 
conductors [l[, extensive studies for the superconduct- 
ing gap structure have been performed because the gap 
structure is closely associated with the pairing mecha- 
nism. Theoretically, fully-gapped s-wave state with op- 
posite signs between different Fermi surfaces (s±-wave) 
has been proposed 0, 0] • The fully-opened gap is sug- 
gested by some experiments, such as penetration depth 
measurements by microwave conductivity 0, @, ARPES 
@, 0], and thermal conductivity Q. However, whether 
the sign-reversal is involved in the multi-gap structure 
is still controversial. While the inelastic neutron mea- 
surements suggests resonance peak in magnetic excita- 
tion spectra x which is observed when the sign 
of the gap takes opposite values on different parts of the 
Fermi surface [^, [T3| , it is pointed out that such a peak 
can be also explained by even s ++ -symmetry, which has 
the same sign of gaps on different Fermi surfaces 11) . 
In addition, several studies on the impurity effect indi- 
cate that the critical temperature T c is robust against the 
introduction of non-magnetic impurities, which is strik- 
ingly different from the suppression of T c predicted in the 
Til - 14 1 . Hence, these results lead us to consider 



s±-wave 

that S4.4. -symmetry should be added to one of the pos- 
sible candidates for the gap symmetry of iron-based su- 
perconductors. Moreover, recent studies suggest that in 
some iron-based superconductors, such as LaFePO [IBj], 
KFe 2 As 2 [l6|, and BaFe 2 (As,P) 2 |17j, the gap is nodal in 
a part of the Fermi surface. In Co-doped BaFe 2 As 2 , it is 
reported that a fully opened gap structure changes to a 
nodal one when the Co-concentration increases from op- 
timal to over-doped region [l8j | . Further studies on the 
gap structure of iron-based superconductors have been 
desired. 

To elucidate the superconducting gap structure, a de- 



tailed study on the effect of defects is very crucial be- 
cause the pair-breaking effects due to scattering centers 
are phase-sensitive. The conventional way to introduce 
impurities is chemical substitutions of constituent ele- 
ments. It is well known that isotropic s-wave supercon- 
ductivity is robust against non-magnetic impurities due 
to Anderson's theorem while superconductivity with a 
sign change in the gap, such as d-wave, is sensitive to 
non-magnetic impurities. However, chemical substitu- 
tions may lead to inhomogeneity in the sample, change 
of carrier density and the Fermi surface topology, which 
can mask the intrinsic impurity effect. Another way to in- 
troduce scattering centers is to create defects by the swift 
particle irradiation. Among them, a light element irra- 
diation, such as proton and a- particle, is very suitable 
for the study of artificially introduced scattering centers, 
since the irradiation can introduce point defects without 
providing inhomogeneity and changing electronic struc- 
ture. So far, for single crystals, only one group reports 
a-particle irradiation experiments F4j. However, the ir- 
radiation produces Kondo-like upturn in the resistivity 
due to spin flip scattering, which can mask the intrin- 
sic non-magnetic scattering effect on the superconducting 
gap of iron-arsenide superconductors. 

In this letter, we address the issue of the supercon- 
ducting gap structure in Ba(Fei_ a; Co :! ;) 2 As 2 (x = 0.045, 
0.075, and 0.113) by detailed study of pair-breaking ef- 
fect introduced by proton irradiation. We find mono- 
tonic increase of the resistivity with proton irradiation. 
The upturn of resistivity at low temperatures is not ob- 
served, which indicates that proton irradiation provides 
non-magnetic scattering centers. The suppression of T c 
is weaker than the expectation for a superconductor with 
sign-reversed gaps on/between the Fermi surfaces. 

Single-crystalline samples of Ba(Fei_ ;r Co :c ) 2 As 2 were 
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FIG. 1: (color online) Temperature dependence of the resis- 
tivity for Ba(Fei_ a; Co :r )2As2 with (a) x = 0.045, (b) 0.075, 
and (c) 0.113. The doses are 0, 0.1, 0.5, 0.8, 1.0, and 1.2xl0 16 
cm" 2 from the lowest curve. The dashed blue lines are fit to 
the data using the equation p — po + AT a with a = 2, 1, and 
1.5 for x = 0.045, 0.075, and 0.113, respectively. 



grown by the FeAs/CoAs self- flux method and their fun- 
damental properties are reported in Ref [[3]]. A mixture 
with a ratio of Ba : FeAs/CoAs =1:5 was placed in 
an alumina crucible. The whole assembly was sealed in a 
large silica tube, and heated up to 1150 °C and kept there 
for 10 h followed by slow cooling down to 800 °C at a rate 
of 5 °C/h, which is slightly different from the synthesis 
reported before [jfj^ . After cleaving, we can obtain shiny 
samples. The typical dimensions of the resulting crystals 
are 4x4x0.1 mm 3 . The average Co concentration in each 
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FIG. 2: (color online) Dose dependence of the normalized 
critical temperature T c /T c o for Ba(Fei_ :E Co: E )2As2 with (a)x 
= 0.045, (b) 0.075, and (c) 0.113. The T c0 is 15.1, 24.8, and 
12.8 K for x = 0.045, 0.075, and 0.113, respectively. Inset: 
Increased residual resistivity by irradiation, Apo = po r — 
Po 



batch was determined by energy dispersive X-ray spec- 
troscopy measurements. The 3 MeV protons, which are 
known to create from one to few tens of displacements 
2l| . were irradiated into the samples at HIMAC-NIRS. 
The irradiation were carried out at 40 K avoiding the 
thermal annealing effect [22j. A total dose is 1.2 x!0 16 



cm . To ensure the uniformity of damage throughout 
the sample, we used samples with thicknesses of 15—30 
/im, which is smaller than the projected range of ^50 
/im obtained from the simulation using the Stopping and 
Range of Ions in Matter-2008 [23| . Resistivity measure- 
ments were performed in-situ after each irradiation by 
standard four-probe configuration. Similar results are 
confirmed in 6 MeV proton irradiation. 

Figure 1 shows the resistivity of Ba(Fei_ a; Co a ;)2As2 
with x = 0.045, 0.075, and 0.113 as a function of temper- 
ature. With increasing dose, T c decreases monotonically 
without significant broadening of the transition width, 
while the resistivity increases monotonically. It should 
be noted that in the a-particle irradiated NdFeAs(0,F) 
Kondo-like resistivity upturn at low temperatures is re- 
ported [3], which is associated with the spin- flip scat- 
tering due to magnetic impurities. In contrast to the 
behavior in NdFeAs(0,F), no upturn is observed in the 
resistivity of proton irradiated Ba(Fei_ :r Co :E )2As2, which 
strongly suggests that defects produced by the irradia- 
tion act as non-magnetic scattering centers. We empha- 
size that only the contribution of non-magnetic scattering 
centers to the pair breaking enable us to investigate the 
intrinsic non-magnetic impurity effect on the order pa- 
rameters of iron-arsenide superconductors. 

Figure 2 shows the dose dependence of the normal- 
ized critical temperature T c /T c q for Ba(Fei_ a; Co a ;)2As2, 
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FIG. 3: (color online) T c as a function of Apo for 
Ba(Fei_ ;c Co :E )2As2 with x = 0.045, 0.075, and 0.113. 



where T c q is the transition temperature before the irradi- 
ation. T c q obtained from the midpoint of resistive transi- 
tion are 15.1, 24.8, and 12.8 K for x = 0.045, 0.075, and 
0.113, respectively. In all the samples, T c /T c q decreases 
linearly with increasing the dose in the present dose 
range. We note that the maximum dose of ~ 1.2 x 10 16 
cm in the present study is one forth of that in the 
a-particle irradiated NdFeAs(0,F). Interestingly, in the 
under- and optimally doped samples, the suppression 
of T c is very small while in the over-doped sample the 
suppression is large down to the half of T c q. We note 
that stronger suppressions of T c in overdoped samples 
are also reported in Zn-doped LaFeAsO0.85F0.15 |24[ and 
LaFeAsOo.85 [Hj]. The increased residual resistivity by 



the irradiation Ap as a function of dose is plotted in 
the inset of Fig. 2. Ap is the difference of the resid- 
ual resistivity between irradiated and unirradiated one, 
namely, Apo — p l J r — pg mrr , which corresponds to the 
density of defects introduced by the proton irradiation. 
We evaluate the residual resistivity po by fitting the data 
using p = pa + AT a , where a is an exponent of temper- 
ature. We fix a as 2, 1, and 1.5 for x = 0.045, 0.075, 
and 0.113, respectively. Apo increases almost linearly 
with dose, which ensures that the proton irradiation in- 
troduces defects systematically. 

Figure 3 shows T c as a function of Apo for 
Ba(Fei_ a; Co a ;)2As2. The suppression of T c due to de- 
fects introduced by the irradiation is almost linear for 
samples with all doping levels. The slope dT c /d(Apo) is 
-0.08, -0.13, and -0.20 K/p,n cm for x = 0.045, 0.075, 
and 0.113, respectively. These values are slightly larger 
than the initial slope of the suppression in the a-particle 
irradiated NdFeAs(0,F), 0.04 K/pi7 cm. 

To discuss the pair-breaking effect due to non-magnetic 
scattering quantitatively, a key parameter is the normal- 
ized scattering rate g — h/2TrkBT C QT. Here, r is the scat- 



0.8 



0.6 



0.4 



0.2 



0.0. 





H + 3 MeV 




Ba(Fe 1 . x Co x ) 2 As 2 




0.045 




-m- 0.075 




0.113 


"•,9c" 




1 , A 





1 2 3 4 5 

g H = f)Apoe/2jtk B T c om"R H 



0.8 



0.6 



0.4 



0.2 



0.0, 





H + 3 MeV 




Ba(Fe 1 . x Co x ) 2 As 2 




-»- 0.045 




-m- 0.075 




0.113 






\9c 




i , ,\ 


"'\ IP) 



1 



0.8 



0.6 



K 
K 



0.4 



0.2 



0.0, 



g 5o ' b = 1 .&ZzApofT cC 





H + 3 MeV 




Ba(Fe 1 . x CoJ 2 As 2 




0.045 




-m- 0.075 




^ 0.113 






1 \ 


(c) 



1 



FIG. 4: (color online) Normalized critical temperature T c /T c q 
as a function of normalized scattering rate (a) g H = 
ftApoe/27rfc s T c0 m*i?H, (b) g 5orb = 1.63^Apo/T c0 , and (c) 
g x = hAp /2irk B T c0 p \ 2 for Ba(Fei_ a: Co a; )2As2 with x = 
0.045, 0.075, and 0.113. Dashed lines are linear extrapo- 
lations. g£ is the critical scattering rate expected in s±- 
scenario. 



tering time including both intra- and inter-band scatter- 
ing contributions. To avoid ambiguity of estimation, we 
present g obtained from three different ways. In order 
to obtain the elastic scattering rate introduced by the 
irradiation, we use the relation t _1 = ne 2 Apo/m* — 
eApo/m*RH, where n is the carrier number and Rh 
is Hall coefficient, m* is the effective mass. Figure 4 
(a) shows T c /T c q as a function of g H = h/2irkBT C QT = 



4 



hAp e/2irk B T m*R H for x = 0.045, 0.075, and 0.113. 
Assuming that electrons are dominant carriers, we use 
the Rh at 300 K obtained from Ref. (2(| and m* ~ 3.5m e 
in the electron pocket obtained from ARPES measure- 
ments for Ba(Fci_ K Co :E )2As2 [27j]. It should be noted 
that the estimated To = m* Rh / epo for unirradated sam- 
ple with x = 0.075 is ~ 0.02 ps, which is consistent 
with the scattering time of ~ 0.05 ps just above T c ob- 
tained by microwave conductivity measurements for K- 
doped BaFe2As2 [H because the residual resistivity for 
unirradated sample with x = 0.075 is about twice as 
large as that of K-doped BaFe2As2- According to the 
s±-scenario with equal gaps of opposite signs on differ- 
ent Fermi surfaces [28], T c obeys the equation,— In t — 
-0(1/2 + g/2t) - ^(1/2), where t = T c /T c0 and ip(x) is 
di-gamma function. This equation indicates that T c van- 
ishes at g = gf < 0.3. To estimate critical values of 
normalized scattering rate for all samples, we linearly 
extrapolated the data for simplicity. The obtained val- 
ues of critical g^ are ~ 6.8, 3.8, and 2.5 for x =0.045, 
0.075, and 0.113, respectively. Even in x = 0.113, where 
T c is most strongly suppressed among them, g^ is much 
larger than expected gf- for s ± -scenario. 

We show another way to estimate g c based on the 
parameter obtained by theoretical calculation. Accord- 
ing to the linear response theory based on five orbital 
model [2^, [3(J, we can obtain the relation Apo[fj,Qcm\ = 
0.098r _1 [K] in Ba(Fei_a; Co^ AS2 with interplane dis- 
tance c = 6.5 A and n = 5.8 — 6.1. Figure 4 (b) 
shows T c /T c q as a function of g 5orb = zTi/2nkBT C QT = 
l^SzApo/TcO: where z = m/m* is the renormalization 
factor. We use z = 1/3.5 obtained by ARPES ^ as- 
suming m ~ m e . Obtained critical values g^ orh by linear 
extrapolation are 6.1, 3.5, and 2.4 for x — 0.045, 0.075, 
and 0.113, respectively, which are again much larger than 
gf for scenario. We note that the critical values g^ orb 
obtained from the parameter based on theoretical calcu- 
lation is very similar to g^ obtained from experimental 
values of Rh and m*. 

To obtain carrier number and effective mass indirectly, 
we use the relation r _1 = ne 2 Ap/m* — Ap/poX 2 , where 
A is the penetration depth, A = ^ p a m* /ne 2 . Figure 4 (c) 
shows T c /T c q as a function of g x = TiApo/2irkBT c opoX 2 . 
Tunnel diode resonator measurements for Al-coated sam- 
ples provides the absolute values of penetration depths 
in Ba(Fei_ x Co 2; )2As2 [Hj]. Above x = 0.045, the ab- 
solute value of penetration depth is almost independent 
of Co doping and is close to 200 nm with small scat- 
tering. Therefore, we use the value of the penetration 
depth A = 200 nm for all samples presented here. Ob- 
tained critical values g x by linear extrapolation are 3.1, 
1.8, and 1.2 for x = 0.045, 0.075, and 0.113, respectively. 
Although these values are roughly half of the previous 
two estimations, they are more than three times larger 
than g^r expected for s ± -scenario. 

Critical scattering rates obtained in the three differ- 



ent estimations in the present study are larger than 
that expected for s± scenario. It should be em- 
phasized that proton irradiation provides only non- 
magnetic scattering centers without changing electronic 
structure of Ba(Fei_ :E Co :E )2As2. Our present results 
definitely indicates that iron-arsenide superconductor 
Ba(Fei_ a: Co 2 ;)2As2 is robust against non-magnetic scat- 
tering. The weak suppression of T c in s±-wave supercon- 
ductors could be understood by the details of scattering 
potential. The pair-breaking effect can be suppressed in 
very weak or negative scattering potential |29j . Strikingly 
suppressed interband scattering introduced by the irra- 
diation, which is rather unlikely in Ba(Fei_ a; Co a ;)2As2 
with both electron and hole pockets having d xz and d yz 
orbital characters, could also explain the weak pair break- 
ing. Further detailed theoretical model to understand the 
non-trivial weak suppression of T c in Ba(Fei_ a; Co a ;)2As2 
should be required. 

Finally we comment on the possibility of change of the 
gap structure with doping level in Ba(Fei_ :c Co :! ;)2As2. 
Thermal conductivity measurements suggest the exis- 
tence of nodes in over-doped sample [18T] . The stronger 
suppression of T c in over-doped sample than under- and 
optimally doped ones may suggest the different _gap struc- 
ture, for instance, nodal s±-wave symmetry [32|, where 
the order parameter has d-wave-like nodes on the elec- 
tron Fermi surface while others are fully open, or s±-wave 
with accidental horizontal nodes [l8| . 

In summary, we present the suppression of T c by 
defects introduced by 3 MeV proton irradiation in 
Ba(Fei_ a; Co a ;)2As2 single crystals with different doping 
levels. We find that T c decreases and residual resistivity 
increases monotonically. No Kondo-like upturn in the low 
temperature resistivity is observed which suggests that 
defects created by the irradiation act as non-magnetic 
scattering centers. The critical scattering rates obtained 
from the three different estimations are much larger than 
that expected in s±-scenario, which may contradict the 
theoretical expectation based on inter-pocket scattering 
due to antiferromagnetic spin fluctuations. 
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